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Abstract 

We study the implications of the recently improved upper limits on the branching 
ratios for the decays B — > /ry, expressed as R(p^/K*'y) = B(B — > fry)/B{B — ► K*~/) < 
0.047. We work out the constraints that the current bound on R(pj/ K*j) implies 
on the parameters of the quark mixing matrix in the standard model (SM). Using 
the present profile of the unitarity triangle, we predict this ratio to be R(p / y/K*^) = 
0.023 ± 0.012. We also work out the correlations involving R(pj/K*'y), the isospin- 
violating ratio A(p^y), and the direct CP-violating asymmetry ^4cp(P7) in B — > 
decays in the SM, in the minimal supersymmetric extension of the SM (MSSM), and 
in an extension of the MSSM involving an additional flavor-changing structure in b — > d 
transitions. 
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1 Introduction 
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Recently, the BABAR collaboration has reported a significant improvement on the upper 
limits of the branching ratios for the decays B°(B°) — > p°7 and B^ — > p ± 7. Averaged over 
the charge conjugated modes, the current 90% C.L. upper limits are [1]: 

(1) 
(2) 
(3) 

They have been combined, using isospin weights for B — > p7 decays and assuming £>(-B° — > 
CJ7) = £>(£>° — > p°7), to yield the improved upper limit 

B(B -> p 7 ) < 1.9 x 10~ 6 . (4) 

The current measurements of the branching ratios for B — > if* 7 decays by BABAR [2], 

B(B° -> iT 7 ) = (4.23 ± 0.40 ± 0.22) x 10" 5 , (5) 
iB(S+ -> iiT* + 7) = (3.83 ± 0.62 ± 0.22) x 10" 5 , (6) 

are then used to set a 90% C.L. upper limit on the ratio of the branching ratios [1] 

This bound is typically a factor 2 away from the SM estimates [3], which we quantify more 
precisely in this letter. In beyond-the-SM scenarios, this bound provides a highly significant 
constraint on the relative strengths of the b — > d'y and b — > 57 transitions. 

The impact of the measurement of R(frf j 'if* 7) on the parameters of the quark mixing 
matrix (henceforth called the Cabibbo-Kobayashi-Maskawa CKM matrix) has been long an- 
ticipated (see, for example, [4]). This quantity measures essentially the CKM matrix element 
ratio |V^| 2 /| V ts \ 2 in the SM. However, one expects significant long-distance contributions in 
R{p r y/K* n l) entering in the decay B — * pj. They are dominated by the annihilation diagrams 
bu — > du-) in the decays B~ — > p~7 [5, 6, 7, 8], which depend on the CKM matrix elements 
V ub V* d . The corresponding annihilation contribution in the decays B° — > p°7 (and its charge 
conjugate) is parametrically suppressed due to the electric charge of the spectator quark in 
B° and the unfavorable color factors. QCD corrections to the decay widths for B — > p7 
also introduce a dependence on V u bV* d in both the charged and neutral 5-meson decays. As 
the relevant CKM matrix element ratio \ u = V u bV* d /VtbV t * d is of 0(1), these modifications 
are important and have to be taken into account in the analysis of R(pj / 'K*j) and other 
observables in B — > p7 decays. 

Recently, the 0(a s ) corrections in the decay widths for B — > V7 (V = K*,p) have been 
calculated in the context of a QCD factorization framework [9], taking into account the 
explicit 0(a s ) and 1/Mb corrections to the penguin amplitudes [3, 10, 11]. Using the theo- 
retical results at hand, we analyze the impact of the current upper limit R(pj / K*j) < 0.047 
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C = 0.76 ±0.10 

A^ K * = -0.113 - i0.043 

A u = -0.0181 + «0.0211 


L U R = -0.095 ± 0.022 
A«* = -0.114 - i0.045 


Vtt = 0.57 

rj tc = 0.47 ±0.04 

i] B = 0.55 

£ s = 1.18 ± 0.04lg 12 


r] cc = 1.38 ± 0.53 

B K = 0.86 ±0.15 

f Bd yB Bd = 235 ± 33 124 M eV 


A = 0.221 ±0.002 

e K = (2.271 ±0.017) 10~ 3 

a ^ Ka = 0.734 ± 0.054 


\V ub /V cb \ = 0.097 ±0.010 
AM Bd = 0.503 ± 0.006 ps" 1 
AM Bs > 14.4 ps" 1 (95% C.L.) 



Table 1: Theoretical parameters and measurements used in B — > p7 observables and in the 
CKM unitarity fits. 



in the context of the SM, where it yields constraints on the CKM parameters p and fj [3], and 
in some popular extensions of the SM, such as the minimal flavor violating minimal super- 
symmetric standard model (MFV-MSSM) [12], and in an Extended-MFV-MSSM scenario 
(EMFV) [13], having a non-CKM flavor-changing structure involving the b — > d transition. 
We also present the correlations involving R(fry / 'if* 7), the isospin-violating ratio A(p7), and 
the direct CP- violating asymmetry A CP (p7) in B — > p7 decays, in the three models just men- 
tioned. Precise measurements of these correlations would provide a strong discrimination 
among the competing models. 



2 Observables 

The effective Hamiltonian for the radiative decays B — > ^7 (equivalently 6 — > c?7 decay) 
can be seen for the SM in [3]. We shall invoke this effective Hamiltonian also for the MFV- 
MSSM and the EMFV cases, which differ from the SM in the Wilson coefficients (WC's), in 
particular in the effective WC's for the magnetic moment operator, Cj (for b — > sj) and Cy 
(for b — > d"f). While this is certainly not the most general operator basis, it is a sufficient 
basis to illustrate the beyond-the-SM effects that may arise in these decays. Restricting 
ourselves to this basis, we first present the 0(a s )-corrected expressions for the observables 
in the B — > p7 decays, worked out in the SM in [3, 10], but now generalized to the case of 
complex Wilson coefficients. 

Here and in the following we will always consider quantities averaged over the charge 
conjugated modes (with the obvious exception of the CP asymmetries). Starting from the 
decay widths T(B+ -> V + -f), T(B~ -> V~ 7 ), T(B° -> V° 7 ) and T(B° -> y° 7 ) (with 
V — p, if*), we construct 

■*(^V7) = ^^ ,l : r( ^ r,l l (8) 
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(9) 



We will define the various observables in terms of these quantities. Note that, up to the 
NLO approximation, this procedure is equivalent to defining two distinct observables for 
the charge conjugate modes and then performing the average. It is preferable to use the 
above definitions since they involve quantities (the CP averaged decay widths) that are 
much easier to measure than the widths of the individual channels, which would require 
tagging the 5-meson. 

The expression for the ratios R(pj / if* 7) is [3] 
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(10) 
(11) 



where ( = £±(0)/£± (0), w hh^(0)(^v5 (0)) being the form factors (at q 2 = 0) in the effective 
heavy quark theory for the decays B — > p(K*)y. Noting that in the SU(3) limit one has 
( = 1, some representative estimates of the SU(3)-breaking and the resulting values of ( are: 
C = 0.76 ± 0.06 from the light-cone QCD sum rules [5]; a theoretically improved estimate 
in the same approach yields [14]: ( = 0.75 ± 0.07; ( = 0.88 ± 0.02 using hybrid QCD sum 
rules [15], and ( = 0.69 ± 10% in the quark model [16]. Except for the hybrid QCD sum 
rules, all other approaches yield a significant SU(3)-breaking in the magnetic moment form 
factors. In the light-cone QCD sum rule approach, this is anticipated due to the appreciable 
differences in the wave functions of the K* and p-mesons. To reflect the current dispersion 
in the theoretical estimates of (, we take its value as ( = 0.76 ± 0.10, given in Table 1. We 
stress that the error (±0.10) is not on ( = 0.76, but rather on the deviation of ( from its 
SU(3) limit, i.e., 1 — ( = 0.24, and amounts to an error of ±42% on the SU(3)-breaking in 
the ratio of form factors in radiative decays. As this is the dominant theoretical error on 
the ratios R ± (prf / K 1 ) and R°(p / y /K*^), it is imperative to reduce it. A lattice-QCD based 
estimate of the form factors, and hence C, is highly desirable. 

The quantity (1 + AR) entails the explicit 0(a s ) corrections, encoded through the func- 



tions A^ K . 

and B ± — > K* ±r y, this can be written after charge conjugated averaging as 
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The definitions of the quantities A^ K \ A^\ A u and L R can be seen in [3]. Their default 
values are summarized in Table 1, where we have also specified the theoretical errors in the 
more sensitive of these parameter L\. The quantity 1 ± AR° is obtained from Eq. (12) in 
the limit L U R = 0. 

The isospin breaking ratio is given by 

r ±( jB - pry) 



A( P7 ) 



2 r°(£ -> pi) 



(13) 
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Figure 1: Unitary triangle fit in the SM and the resulting 95% C.L. contour in the p - rj 
plane. The impact of the R{p r y/K* , -y) < 0.047 constraint is also shown. 
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and the CP asymmetry in the charged modes is 



Acp{pl) 
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The CP asymmetry in the neutral modes 
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(14) 



(15) 
(16) 



(17) 



is obtained from Eq. (16) in the limit L U R = 0. 

Note, that as in the EMFV model there are additional contributions to the effective 
Wilson coefficients C7, entering through the mass insertion parameters (for the b — > d 7 
case) and 5% 3 (for the b — > s 7 case), which are in general different, we have introduced two 
different magnetic moment WC's for the d and s sectors. For Cj = Cj = Cj M we reproduce 
the formulae presented in Ref. [3] with the only exception of Ai?(p 7 ); in this case, the factor 
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Figure 2: Extremal values of R(pj/K*j) that are compatible with the SM unitarity triangle 
analysis. 

in the parentheses in the first line of Eq. (12) is missing in Ref. [3]. This, however, has 
only a small numerical effect, as can be judged from the values AR = 0.055 ± 0.130 and 
AR° = 0.015 ± 0.110 that we have obtained here and which are in quite good agreement 
with the corresponding values given in Ref. [3]. 



3 Impact on the unitarity triangle 



In this section, we present an updated analysis of the constraints in the (p, fj) plane from 
the unitarity of the CKM matrix, including the measurements of the CP asymmetry a^K s m 
the decays B°/B° — > J/ipK s (and related modes), and show the impact of the upper limit 
R{frrfK*i) < 0.047 [1]. 

The SM expressions for ex (CP- violating parameter in K decays), AM Bd (B®-B® mass 
difference), AMg s (B®-B® mass difference) and a^K 8 are fairly standard and can be found, 
for instance, in Ref. [17], where also references to the various theoretical input parameters 
which have not changed since then can be found. Note that for the hadronic parameters 

f Bd \J B Bd and £ s , we use the recent lattice estimates [18] which into account uncertainties 
induced by the so-called chiral logarithms [19]. These errors are highly asymmetric and, 
once taken into account, reduce sizeably the impact of the AM B J AM Bd lower bound on 
the unitarity triangle analysis. The experimental inputs for the quantities A and ex are 
taken from the Particle Data Group [20]. The measurement of the CP asymmetry a^x a i n 
the decays B°/B° — > J/ipK s (and related modes) is now dominated by the BABAR [21] 
and BELLE [22] collaborations; taking into account the earlier measurements yield the 
current world average a^K s — 0.734 ± 0.054 [23]. The indicated value of the mass difference 
AM Bd = 0.503 ± 0.006 ps~ x is the current world average [24] and the 95% C.L. lower bound 
AM Bs > 14.4 ps _1 has been recently updated this summer [25]. The values of the theoretical 
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parameters and experimental measurements that we use are summarized in Table 1. 

The SM fit of the unitarity triangle is presented in Fig. 1, where we show explicitly what 
happens to the allowed regions once the errors associated with the chiral logs are taken into 
account. The 95% CL. contour is drawn taking into account chiral logarithm uncertainties. 
The fitted values for the Wolfenstein parameters, the angles a and (3, AM Bs and the CKM 
ratio iV^/Kifcl are given below where we also show the resulting values that we obtain without 
including the chiral logarithms uncertainties: 



P 
fj 
a 

P 
1 

AM Bs 



X-logs 
0.22 ±0.07 
0.34 ±0.04 
(98 ± 10)° 
(24.2 ± 1.8)° 
(60 ± 10)° 

(19.61ft) ps" 
1.75 ±0.15 



no x~l°g s 
0.25 ±0.07 
0.34 ±0.04 
(101 ± 10)° 
(25.0 ± 1.9)° 
(56 ± 10)° 
(21.0if;t) PS" 
1.61 ±0.14 



The main effect of the chiral logs is that they decrease the central value of p by about half 



a sigma with fj remaining practically unchanged. The 
increased value of the CKM matrix element ratio 
about 1 sigma. 



Vtd 



V,, 



argest impact of this shift is in the 
, whose central value moves up by 



As the bound from the current upper limit on i?(p7 / 'if* 7) is not yet competitive to the 
ones from either the measurement of AM Bd , or the current bound on AM Bs , we use the 
allowed p — fj region in order to work out the SM predictions for the observables in the 
radiative 5-decays described above. Taking into account these errors and the uncertainties 
on the theoretical parameters presented in Table 1, we find the following SM expectations 
for the radiative decays: 



R ± (p 1 /K* 1 ) 
R Q {p 1 /K* 1 ) 
A(p 7 ) 

^•p(pl) 
A° CP (P1) 



0.023 ±0.012 , 
0.011 ±0.006 , 
04 +0 - 14 
10+ - 03 

u - iu -0.02 ) 

0.06 ±0.02 . 



(18) 
(19) 
(20) 
(21) 
(22) 



It is interesting to work out the extremal values of R(pj/ 'if* 7) compatible with the SM UT- 
analysis. This is geometrically shown in Fig. 2 where we draw the bands corresponding to the 
values 0.038 ±0.009 and 0.013 ±0.003 (the errors are essentially driven by the uncertainty on 
(). The meaning of this figure is as follows: any measurement of R(pry j 'if* 7), whose central 
value lies in the range (0.013, 0.038) would be compatible with the SM, irrespective of the 
size of the experimental error. The error induced by the imprecise determination of the 
isospin breaking parameter ( limits currently the possibility of having a very sharp impact 
from R(pj/K* r y) on the UT analysis. 
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Figure 3: Correlation between R(p r y/K* r y) and A(p7) in the SM and in MFV and EMFV 
models. The light-shaded regions are obtained varying p, fj, the supersymmetric parameters 
(for the MFV and EMFV models) and using the central values of all the hadronic quantities. 
The darker regions show the effect of ±la variation of the hadronic parameters. 
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Figure 4: Correlation between R^pj / K*^) and A CP (fr/). See the caption in Fig. 3 for further 
details. 

4 Analysis in supersymmetry 

We focus on two variants of the MSSM called in the literature as MFV [12] and Extended- 
MFV [13] models. In MFV models all the flavor changing sources, other than the CKM 
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matrix, are neglected and the remaining parameters (that are assumed to be real) are the 
common mass of the heavy squarks other than the lightest stop (Mq) , the mass of the lightest 
stop (Mi), the stop mixing angle (9-), the ratio of the vacuum expectation values of the two 
Higgs bosons (tan (3 S ), the two parameters of the chargino mass matrix (// and M 2 ) and 
the charged Higgs mass (M H ±). In this class of models there are essentially no additional 
contributions (on top of the SM ones) to cl$k s and AMbJ AMg^, while the impact on 
€k, AMe d and AMg s is described by a single parameter /, whose value depends on the 
parameters of the supersymmetric models [17]. 

In EMFV models, there is an additional parameter 

M 2 r v M 

M iMq \V td \ ■ [ 6) 

With the inclusion of this new parameter, the description of the UT-related observables 
needs one more complex parameter, g = ga + igi [13]. A signature of these models is 
the presence of a new phase in the B® — B^ mixing amplitude. Using the parametrization 
Mf 2 = rje 2i9d Mf^, we get rj = |1 + / + g\ and 9 d = l/2arg(l + f + g). This implies 
supersymmetric contributions to the CP asymmetry a^K s , which we quantify below. 

We analyze the phenomenology of the MFV and EMFV models by means of scatter plots 
over the supersymmetric parameter space. In the MFV case, we scan over the following 
ranges (M q is set to 1 TeV, likewise M- g is 0(1 TeV)): 

(24) 
(25) 
(26) 
(27) 
(28) 

In the EMFV case, we limit the range of the stop mixing angle to Oi = [—0.3 -j- 0.3] (see 
discussion in Ref. [13]) and add the scan over f| = [0^-1] and axgS^^ = [—n -v- n]. 
We scan also over p and fj and require that each point satisfies the bounds that come from 
direct searches, from the B — > branching ratio, and from the UT related observables 
summarized in Table 1. The surviving regions are presented in Figs. 3-5. In each figure, 
the light shaded regions are obtained using the central values of the input parameters given 
in Table 1 while the dark shaded ones result from the inclusion of their one sigma errors. 
Note that in the two figures showing the correlations between Aq-p(p^) and R(frf/K*^), and 
^cv(Pl) an d R{P1 1 respectively, the CP asymmetries tend to increase as expected in 
the limit of small branching ratios. In the MFV case, there are two distinct regions that 
correspond to the negative (SM-like) and positive C| cases. For Cj < 0, the allowed regions 
in MFV almost coincide with the SM ones and we do not draw them. For Cj > 0, the allowed 
regions are different and, in general, a change of sign of both the CP-asymmetries (compared 
to the SM) is expected. We note that the latter scenario needs very large SUSY contributions 
to Cj, arising from the chargino-stop diagrams, and for fixed values of tan /3s it is possible 
to set an upper limit on the mass of the lightest stop squark. In Fig. 6, we show the points 
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Figure 5: Correlation between R{p r )/K* r y) and A° CP (p r )). See the caption in Fig. 3 for further 



that survive the B — > X s "y constraint with a positive C| in the MFV scenario. We have also 
imposed the additional constraint coming from the upper limit B(B S — > < 2.6 x 10 -6 

at 90% C.L. [27], and find that the allowed region is essentially unaffected. Note that in the 
> scenario the mass of the lightest stop has an upper bound of 500 GeV for tan j3s < 50. 

What concerns the allowed values of the phase 9 d , comparing the SM allowed range from 
the UT fit sin 2/3 = 0.76 ± 0.06, implying (3 = 25° ± 3°, with the current experimental value 
a %l>K a = 0.734 ± 0.054, yields 9 d G (—5°, 8°). The other solution for a^ Ks shown in Fig. 1 
yields 9 d G (33°, 46°). The additional contributions in Mf 2 impact on the dilepton charge 
asymmetry [26] 



where £ ++ {£ ) are the numbers of £ + £ + (£~£~) observed in the decay of a BB pair, and Ar^ 
is the difference in the decay widhts of the two mass eigenstate. We have computed for each 
point the value of the dilepton asymmetry Am and found that the allowed range in the EMFV 
model is Au G [—0.1,0.7] x 10 -2 . This expectation has to be compared with the current 
experimental bound, A^ p = (0.46 ± 1.18 ± 1.43) x 10~ 2 [28]. We see that the experimental 
precision has to be improved by one order of magnitude in order to test the EMFV models. 
Note that (AT d /AM d ) SM = (1-3 ± 0.2) x 10" 2 , yielding typically A u = O(10" 3 ) in the SM. 
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Figure 6: Correlation between tan (5s and Mi (in GeV) in the MFV-SUSY model for the 
Cf > scenario. 

5 Summary 

We have presented here an analysis of the ratio R(py/K*^), involving the decays B — > p7 
and B — > K*j, the isospin-violating asymmetry in B — > p7 decays A(p7), and direct CP 
asymmetries A^p(p7) and Ap P (p7) in the charged and neutral 5-meson decays in the SM 
and two variants of supersymmetric theories. They illustrate the current and impending 
interest in the radiative decays B — > p7, which will provide powerful constraints on the 
CKM parameters and allow to search for physics beyond-the-SM. 
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